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ABSTRACT 

We assume that the hehum-I hnes emitted by the massive binary system rj 
Carinae are formed in the acceleration zone of the less-massive secondary star. 
We calculate the Doppler shift of the lines as a function of orbital phase and 
of several parameters of the binary system. We find that a good fit is obtained 
if the helium lines are formed in the region where the secondary wind speed is 
Vzone = 430 km s~^. The acceptable binary eccentricity is in the range 0.90 < 
e < 0.95, and the inclination angle (the angle between a line perpendicular to the 
orbital plane and the line of sight) is in the range 40° ^ i ^ 55°. Lower values 
of e require higher values of i, and vice versa. The binary system is oriented 
such that the secondary star is in our direction (closer to us) during periastron 
passage. The orbital motion can account in part to the Doppler shift of the peak 
in X-ray emission. 

Subject headings: (stars:) binaries: general— stars: mass loss— stars: winds, 
outflows— stars: individual {rj Carinae) 



1. INTRODUCTION 

If it was not for its binary companion (the secondary rj Car B star), the primary massive 
star {i] Car A) would be a 'boring' star. It is the binary companion that most likely shaped 
the Homunculus (Soker 2003; 2007b)— the bipolar nebula around rj Car that was ejected 
during the twenty years long Great Eruption (Davidson & Humphreys 1997). The orbital 
motion is also behind the 5.54 yr year periodic variation in the light curve of r] Car, from the 
radio (Duncan & White 2003) to the X-ray band (Corcoran 2005), and in all wavelengths in 
between (e.g., Damineh 1996; Damineh et al. 1998, 2000; Whitelock et al. 2004; Smith et 
al. 2004; Abraham et al. 2005; Hillier et al. 2006; Nielsen et al. 2007). 
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For some emission bands, and emission and absorption lines, there is a consensus as to 
their origin, for other there is a dispute. It is agreed that the hard X-ray emission results 
from the hot postshock secondary stellar wind (Corcoran et al. 2001; Pittard & Corcoran 
2002; Akashi et al. 2006; Hamaguchi et al. 2007), that most of the radio emission comes 
from extended (~ W^AU) circumbinary ionized gas (Cox ct al. 1995; Duncan & White 
2003; Kashi & Soker 2007), and the source of the hard UV radiation is rj Car B (Verner et 
al. 2005; Duncan & White 2003; Soker 2007a). The location of the plasma emitting the 
He II A4686A line is in dispute. Steiner & Damineli (2004) and Martin et al. (2006) assumed 
that the He II emission originates from photoionized regions near the X-ray shock fronts of 
the two winds. Soker & Behar (2006), on the other hand, postulate that this hue originates 
in the acceleration zone of the secondary wind. 

In the present paper we dispute the recent suggestion made by Nielsen et al. (2007) 
that the Hel lines originate in the primary wind. This is described in section 2. In section 
3 we describe the blueshift resulting from the orbital motion. In section 4 we analyze the 
velocity of the lines and try to constrain the orientation angle of the semimajor axis and 
the inclination of the orbital plane. The orientation angle of the semi-major axis of rj Car is 
another point of dispute. We summarize in section 5. 

2. THE HELIUM-I VISIBLE LINES 

The He I lines (A5877, 6680, 7067, 7283) are formed by transitions between excited 
levels in neutral helium. The main properties of these lines in rj Car are as follows (Nielsen 
et al. 2007). The lines are observed both in absorption and emission, and have P-Cygni 
type profiles at many orbital phases of r] car. The velocity measured at the minimum 
of the line profile (maximum absorption) relative to the system velocity along the line of 
sight is blueshifted at all phases. After periastron passage the velocity is ~ —300 km s~^. 
The magnitude of the velocity then increases, reaching ~ —500 km s~^ at phase ~ —0.05 
(~ 100 days before periastron). Very close to periastron passage there is a sharp rise in the 
magnitude of the velocity, probably due to a second wind component (Nielsen et al. 2007), 
to ~ —600 km s~^. During the spectroscopic event the absorption line velocity rapidly shifts 
from ~ —600 km s~^ to ~ —300 km s~^. 

The equivalent widths of the lines in emission are more or less constant along most 
of the orbit (beside some fluctuations). At phase ~ —0.1 they start to drop, reaching a 
minimum at phase ~ 0, where they stay for the spectroscopic event of ~ 10 weeks. The 
equivalent widths in absorption are low during most of the orbit. At phase ~ —0.1 they 
start to rise reaching a maximum just before phase zero. During the spectroscopic event they 
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sharply drop, partially recovering after the event (for more detail see Nielsen et al. 2007). 
The general behavior of the Hel lines in absorption as reported by Nielsen et al. (2007) 
is qualitatively similar, in both velocity and equivalent width, to the behavior of the He II 
A4686 emission line as reported by Steiner & Damineli (2004; see also Martin et al. 2006). 

Nielsen et al. (2007) proposed that the Hel emission and absorption originate in the 
primary wind, where the helium is ionized (or excited) by the secondary stellar radiation. 
Luminous blue variable stars similar in properties to t] Car A are known to have Hel lines 
with P-Cygni type profiles, e.g., P-Cygni itself (Stahl et al. 1993; Najarro et al. 1997), 
without any external ionizing source. Therefore, it is possible in principle that the source 
of the Hel lines is indeed the primary wind. However, there are some difficulties with this 
possibihty: (1) What cause the monotonic variation in the blueshift along the orbit? (2) Why 
do the equivalent widths of absorption and emission drop for several weeks near periastron 
passage? (3) The most severe problem is the behavior of the He I emission. Not only the 
absorption, put the entire emission part changes its Doppler shift with phase. The width 
of the emission peak does not change much. This is expected if the line is formed by the 
orbiting secondary star, as we propose. To the contrary, if the line originates in the primary 
wind, then as regions further out become the source of the line, with their higher velocity, 
so is the emission peak suppose to be wider. This is not observed (fig. 2 of Nielsen et al. 
2007). Considering these difficulties, we follow Soker Sz Behar (2006) who attributed the 
He II A4686 fine to the wind from the t] Car B, and propose the same for the Hel lines. 
Namely, we argue that the He I lines originate in the acceleration zone of the secondary 
wind. 

WR and O stars are known to have Hell and Hel lines (e.g., Crowther & Bohannan 
1997). The three stars studied by Crowther & Bohannan (1997) are somewhat cooler than 
rj Car B, but the main factor is the mass loss rate (Crowther & Bohannan 1997), which 
is similar to that of rj Car B. These three stars all show the He II A4686 line in emission, 
and the He I A5876 line in emission and absorption with equivalent widths of ~ 0.25 — 4 A. 
These stars have He 1 A6678 emission, and one has this line in absorption as well (Crowther 
& Bohannan 1997). The minimum (maximum absorption) in the P-Cygni line profiles of the 
Hel fines of O stars occurs typically at a velocity of ~ —0.1 to —0.3 times the wind's terminal 
speed (Prinja et al. 1996, 2001, 2006). For the r] Car B this range is ~ 300 - 900 km s-^ in 
accord with the observations of Nielsen et al. (2007). We conclude that the Hel fine profiles 
of O stars are similar enough to the line profiles observed from 77 Car that we may attribute 
these lines to the wind from the secondary star. 

The fast temporal variability in the intensity and speed of these fines from 0-stars (e.g. 
Prinja et al. 2006) is evidence of their sensitivity to the wind properties. Therefore, small 
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changes in the wind properties might resuh in large changes in the observed hnes profiles. 
Soker & Behar (2006) attributed the large changes in the properties of the He II A4686 
emission line to slight changes in the density of the secondary stellar wind in its acceleration 
zone. The increase in density might result from the influence of the enhanced X-ray emission 
arising from the shocked secondary wind further downstream, and/or from accreting mass 
from the primary stellar wind. 

A monotonic change after periastron passage might result from the accretion event. Al- 
though not much mass is accreted during the ~ 10 week long accretion phase, it is enough to 
change the structure of the outer envelope of the secondary star. The secondary itself is still 
recovering from the 19th century Great Eruption, when it accreted more than lOM© (Soker 
2007b). Therefore, it is possible that the secondary requires several years to reestablish its 
structure, and before it managed to do so it encounters the next accretion phase. 

Soker & Behar (2006) speculated that the high velocity component appearing very close 
to periastron passage results from accreted primary wind material by the secondary star. At 
the beginning of the accretion phase the accreted gas has a relatively high specific angular 
momentum, and might coUimate fast segments of the secondary wind. This might explain 
the new wind component mentioned by Nielsen et al. (2007) to appear near periastron. 

Despite the above description, it will be much easier to account for the velocity change 
along the orbit if we observe rj Car such that the secondary is closer to us near periastron, 
namely, 7 ~ (see Figure [T]) . This implies that the secondary moves away from us after 
periastron and toward us before periastron. This is opposite to the orbit orientation sug- 
gested by Nielsen et al. (2007). In both cases the orbital plane is tilted by ~ 45° to the line 
of sight, e.g., Davidson et al. (2001) give the value i ~ 41°, where i is the angle between a 
line perpendicular to the orbital plane and the line of sight. We therefore turn to study the 
blueshift expected from lines emitted or absorbed in the acceleration zone of the secondary 
wind. 



3. LINES VELOCITY 

The 7] Car binary parameters used by us are as in the previous papers in this series 
(Soker 2005; Akashi et al. 2006; Soker & Behar 2006; Soker 2007a), and are compiled by 
using results from several different papers (e.g., Ishibashi et al. 1999; Damineli et al. 2000; 
Corcoran et al. 2001, 2004; Hillier et al. 2001; Pittard & Corcoran 2002; Smith et al. 2004; 
Verner et al. 2005). The assumed stellar masses are Mi = I2OM0, M2 = 30Mq, and orbital 
period 2024 days, hence the semi-major axis is a = 16.64 AU. The winds properties are not 
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directly relevant to us, but the secondary wind mass loss rate and velocity are similar to 
other stars having the Hell and Hel lines (see section 2); these are M2 — 10~^Mq yr~^ and 
V2 — 3000 km s~^, respectively. The binary system configuration is presented in Figure [H 
while in Figure [2] we present the variation of the relative speed between the two stars Vo^h, of 
the orbital separation r, and of the orbital angle 6, with phase (equals zero at periastron), 
for e = 0.9. The angle 9 is the relative direction of the secondary to the primary as measured 
from periastron. 

The blueshifted line velocity is given by 

^line ^zone (1) 

where v^one is the velocity of the wind's zone where the line forms relative to the secondary 
star, and Vm is the modulation resulting from the orbital motion. For a given orbital period 
and orientation angle 7, the blueshifted line velocity dependance on the binary masses, 
inclination angle i, and a large eccentricity e is given by 

t^..oc^^^^-^|-^(l-e)-/^inz. (2) 

the dependance on the masses and inclination angle is accurate, while the dependance on 
eccentricity is a good approximation for large eccentricity. In the lower row of Figure [2] we 
present the line of sight velocity of a spectral line formed where v^one = 430 km s~^, for 
i = 45°, e = 0.9, and 7 = 0. 



4. FITTING THE LINE VELOCITY CURVE 

In this section we change the orientation angle 7 (see Figure 1), the inclination angle 
i, and the eccentricity e to fit the results of Nielsen et al. (2007). In the figures below the 
fitted line velocity from Nielsen et al. (2007) is drawn as a solid thick line. We found that in 
all case v^one = 430 km s~^ (see eq. [1]) gives the best results, and we will use this value in all 
cases. In principle we could play also with the masses of the two stars, hence with the orbital 
semimajor axis. However, because the line speed goes as ~ (eq. [2]), for a small change 
in the line speed a large change in the primary mass is required, mainly higher mass than 
Ml = I2OM0 is required, while most papers cite this or lower masses (e.g., Corcoran 2007). 
Lower secondary mass than we use is also problematic (Verner et al. 2005). Therefore, we 
don't change the masses of the stellar components. 

We fitted the data points from Nielsen et al. (2007), rather than the curve they draw 
through their points. However, for clarity, in the figures below we draw their curve instead 
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secondary wind x-ray 



Fig. 1. — The orbital plane and a cut through the primary wind shocked region; this region 
extends beyond the size of the figure, as drawn schematically at phase t = —0.06. The phases 
are marked at four points along the orbit. The elliptical orbit, the primary stellar size, the 
distance of the shocked primary wind region and its width, are drawn approximately to scale 
for e = 0.9 (but not the secondary star or the secondary wind shock region). The thick 
section of the orbit is approximately where accretion phase occurs according to the model 
(Soker 2005; Akashi ct al. 2006). At phase t — 0.5 we mark the sources of the X-ray emission 
(the shocked secondary wind), of the ionizing radiation (the secondary photosphere), and 
of the Hel and Hell lines (the acceleration zone of the secondary stellar wind). The orbital 
plane is tilted by i ~ 45° to the plane of the sky. 
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Phase Phase 

Fig. 2. — Several physical variables as a function of orbital phase (phase zero is at periastron), 
for e = 0.9. The left column covers the entire orbit, while the column on the right covers the 
time just prior to and after periastron passage. Top row: Orbital separation (in AU) and 
relative orbital speed of the two stars (in 10 km s~^ on the left and 100 km s~^ on the right). 
The angle 9 is the relative direction of the two stars as measured from periastron (scale on 
the right in degrees). Bottom row: The variation of the blueshifted line speed with phase 
(in km s^^. The line is formed in the acceleration zone where the secondary wind speed 
relative to the secondary star is v^one — 430 km s~^, the orientation angle is 7 = 0, and the 
inclination is i = 45°. 
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of their data. For that, the fit does not appear perfect. In addition, we will not give high 
weight to an exact fitting very close to periastron passage. The reason is that very close to 
periastron the properties in the acceleration zone of the secondary wind can change because 
of the X-ray emission and accretion (Soker & Behar 2006). Soker & Behar (2006) argued 
that such property variations can cause the He II A4686 line to be formed in a region with a 
higher velocity near periastron. For that, in all graphs below there is no perfect fitting very 
close to phase zero. 

We start by taking i = 41° (Davidson et al. 2001), e = 0.9 as was used in previous 
papers in this series, and 7 = 0. As evident from the thin solid line in Figure [3l these 
parameters give line velocity with amplitude lower by ~ 30 km s~^ than the line velocity 
found of Nielsen et al. (2007). When comparing the observations (Nielsen et al. 2007) over 
the entire period, we find that for these parameters our modelled line velocity is lower by 
~ 20 km s^^ than required. By varying the orientation of the semimajor axis to 7 ~ 10° 
or 7 ~ —10° we can make a better fitting either before or after the spectroscopic event, 
respectively, but not in both epochs. 

Next we try to increase the inclination angle i, while keeping the other parameters 
unchanged. The dependance of the velocity amplitude on i is according to equation ([2]). We 
find the best fit to be for i ~ 53°. This is somewhat larger than the value considered in some 
previous papers (e.g. Davidson et al. 2001), but not in others (e.g., Corcoran 2007). The 
line velocity for this case is drawn in Figure HI 

Another option to get a better fit is to increase the eccentricity e. Keeping all parameters 
unchanged, namely, i = 41° etc., we find the best fitting for e ~ 0.93, as expected from eq 2. 
The result for this case is plotted in Figure [51 We also try an intermediate case, with values 
of e = 0.92 and i = 45°. The result for this case is plotted in Figured 

In all cases we find that even a small variation as 20° from 7 = decreases the fitting. 
We conclude that the semimajor axis orientation angle is —10° ^ 7 ^ 10°. 

5. DISCUSSION AND SUMMARY 

We postulated that the Hel lines are formed in the acceleration zone of the secondary 
wind, rather than in the primary wind as suggested by Nielsen et al. (2007). Supporting 
arguments for this assumption are given in Section 2. Under this postulate, the Doppler 
shift of the Hel lines and its variation with the orbital phase is caused by the orbital motion 
of the secondary around the center of mass of the binary system. We found that using 
binary parameters from the literature we can nicely fit the Doppler shift. In section 3 we 
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Fig. 3. — The Hel lines velocity fitted by Nielsen et al. (2007; solid thick line) is compared 
with the model discussed here for an inclination angle (between the orbital axes, which 
is perpendicular to the orbital plane, and line of sight) of i = 41° and an eccentricity of 
e = 0.9. Calculations for three orientation angles (see Figure 1) are presented with thin 
lines: 7 = —20° (dotted); 7 = (sohd); 7 = 20° (dashed). The left side covers the entire 
orbit, while that on the right covers the time just prior to and after periastron passage. 




Fig. 4.— Like Fig. [3l but for i = 53° and e = 0.9. 




Fig. 6.— Like Fig. El but for i = 45° and e = 0.92 
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described the binary system and the calculation of the Doppler shift. In section 4 we fitted 
the calculated Doppler shift to the lines velocity observed by Nielsen et al. (2007) by varying 
the inclination angle i, the eccentricity e, and the orientation angle of the semimajor axis 
7. In all cases we found that the best fit is obtained if the Hel lines are formed in the 
acceleration zone of the secondary wind where the wind speed is v^one = 430 km s~^ (see eq. 

CD. 

In all cases the best fit is for an orientation angle 7 = 0. Namely, the secondary is 
closer to us at periastron passage with —10° < 7 < 10°. Falceta-Goncalves et al. (2005) 
already argued for this orientation. However, they based their claim on a model for the 
X-ray light curve which we disagree with (Akashi et al. 2006). Others (e.g., Nielsen et al. 
2007; Corcoran 2007) prefer the orientation to be opposite to our finding, namely, 7 = 180°. 
Smith et al. (2004) used e = 0.8, and argued for a perpendicular orientation with 7 = 90°. 
The orientation angle is one of the binary parameters in high priority for the community to 
agree upon. 

According to our results, and considering the several sources of uncertainties (both in 
model and observations of the lines velocities), the eccentricity is in the range 0.90 ^ e < 
0.95, and the inclination angle in the range 40° ^ i ^ 55°. Lower values of e require higher 
values of i. Because the inclination is constraint by other observations (e.g., Davidson et la. 
2001) to z ~ 41 — 45°, we prefer the range 0.92 ^ e < 0.95 This eccentricity is higher than 
the value of e = 0.9 used in the papers dealing with the accretion model near periastron 
(Soker 2005; Akashi et al. 2006), and implies closer approach of the two stars by > 20%. 
This substantially increase the likelihood of the accretion model. 

We have the following answers to comments raised by some colleagues against our as- 
sumption that the He I lines originate in the secondary wind. (1) Some hydrogen lines with 
similar profiles (Nielsen et al. 2007) originate in the primary wind. Answer: We note that 
the P-Cygni hydrogen line profile, e.g., H I A4103, does not change as the He I lines do 
(Nielsen et al. 2007). (2) The coincidence that the secondary wind velocity where the lines 
are formed, fzone = —430 km s~^, is practically the same as the primary wind terminal 
velocity. Answer: Although this is a coincidence, the same can be said on the model where 
the He I lines originate in the primary wind. How come the region where the lines are 
formed in the primary wind changes its velocity exactly as the secondary speed around the 
center of mass does? (3) Other emission lines expected from the secondary are not observed, 
e.g., the He II A4686 line. Answer: the secondary wind must be highly clumpy, with dense 
and somewhat cooler blobs emitting the He I lines. The sensitivity of lines intensities to 
the exact wind properties can account for the non-presence of other expected lines. (4) 77 
Car B is much fainter than r] Car A, and it would be difficult to see the emission from B. 
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Answer: The secondary is hot and its luminosity is ~ 20% that of the primary. This seems 
enough to account for the He I hues from the secondary to be detected. We note that in 
the model of Nielsen et al. (2007), the secondary will ionize a small section of the primary 
wind, and it is hard to explain the strong emission. (5) The secondary radiation will ionize 
the primary wind, and we expect He 1 lines from the primary wind. Answer: We agree with 
that. However, we expect the contribution from the primary wind to be small, as stated in 
point (4) above. 

What we have shown here is that if the lines originate in the acceleration zone of the 
secondary wind, then their Dopplcr shift can be easily explained. This is not the case if 
the lines originate in the primary wind. In any case, we stress again that He 1 lines are 
expected from r] Car B. PuUerton et al. (1996) listed several stars with effective temperature 
of ~ 36, 000 — 40, 000 K that have the He I 5876 line in absorption. The line velocity is 
typically much lower than the wind terminal speed. These lines are formed in the wind; for 
example, Rauw et al. (2001) mentioned that the Hel lines are formed in the wind of HD 
192639, which has an effective temperature similar to that of r) Car B (FuUerton et al. 1996). 

Another possibility is that the He I lines originate near the stagnation point, namely, 
in the post shock primary wind. The stagnation point speed around the center of mass is 
~ 70% that of r] Car B. This possibility requires a higher eccentricity, e ~ 0.95 for i = 53°, 
or a higher inclination, i ~ 70° for e = 0.93, or some combined increased in their values. 

If the binary parameters we found here are correct, then they have implications for 
Doppler shift in other wavelengths as well. Consider X-ray emission. Most of the hard X-ray 

emission comes from the shocked secondary wind near the stagnation point of the colliding 
winds (Pittard et al. 1998; Corcoran et al. 2001; Pittard & Corcoran 2002; Akashi et al. 
2006; Hamaguchi et al. 2007). This region is close to the secondary star and moves with 
it. Therefore, a Doppler shift is expected in the X-ray band. Indeed, very near periastron 
passage, phase < 0.03, the peak in the X-ray spectrum is blueshifted by ~ 150 — 250 km s^^ 
(Behar et al. 2006). No blue shift is observed about a year before periastron. We therefore 
suggest that this shift in the peak is due mainly to the orbital motion of the secondary. Some 
contribution can come from a second wind component appearing near periastron (Behar et 
al. 2006). Corcoran (2007) took 7 = 180° and proposed that the changes in the Doppler 
shift of the X-ray lines with orbital phase result from the change in the orientation of the 
colliding winds region. In principle such a model can work when we observe the system in 
the orbital plane, namely i = 90°. However, if we take the inclination into account, then even 
during apastron passage a blueshift component of Va — ~ Vp cos i is expected, where Vp is the 
(blueshifted) velocity causing the blueshift the model supposes to explain near periastron. 
Taking i < 55° for 77 Car, we find that Va > 0.5vp. Such a blue shift is not observed near 



- 13 - 



periastron (Behar et al. 2006; Corcoran 2007). 
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